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Abstract 

We introduce a dynamical description based on a probability density (f){a, x, y, t) 
of the vorticity a in two-dimensional viscous flows such that the average vor- 
ticity evolves according to the Navier-Stokes equations. A time-dependent 
mixing index is defined and the class of probability densities that maximizes 
this index is studied. The time dependence of the Lagrange multipliers can 
be chosen in such a way that the masses m{a,t) := Jdxdy (j){a,x,y,t) asso- 
ciated with each vorticity value a are conserved. When the masses m{a, t) 
are conserved then 1) the mixing index satisfies an H-theorem and 2) the 
mixing index is the time- dependent analogue of the entropy employed in 
the statistical mechanical theory of inviscid 2D flows [Miller, Weichman & 
Cross, Phys. Rev. A 45 (1992); Robert & Sommeria, Phys. Rev. Lett. 
69, 2776 (1992)]. Within this framework we also show how to reconstruct 
the probability density of the quasi-stationary coherent structures from the 
experimentally determined vorticity-stream function relations and we pro- 
vide a connection between this probability density and an appropriate initial 
distribution. 



1 Introduction 

When studying the dynamics of two-dimensional fluid motion characterized 
by a vorticity field u!{x,y,t) it can be useful to turn to a statistical descrip- 
tion with probabihty distributions 0((T, x, y, t) for the microscopic vorticity a 
such the average value of a over these distributions is equal to uj{x,y,t). In 
particular, this has been done in the description of the quasi-stationary states 
(QSS), i.e., the coherent structures which are often reached in (numerical) 
experiments after a fast mixing process has taken place [^[TTlfTKt fT ^ I^. At 
high Reynolds' numbers, the vorticity fields ujs{x,y) of these QSS's satisfy 
u — relations to a good approximation, i.e., ujs{x,y) — VL^ip^x^y)) where 
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il){x,y) is the corresponding stream-function. In other words, the QSS's are 
approximate stationary solutions of the Euler equation. 

More specifically, in the early 1990 's Miller [16] and Robert [211 [22] together 
with their coworkers [T71 IMl [23] presented a statistical mechanical theory of 
steady flows in inviscid, two-dimensional fluids, an approach that can be 
traced back to Linden-Bell's work of 1967 [13] • Some outstanding aspects of 
this non-dissipative system are: 1) an infinite number of conserved quantities, 
the masses associated with each microscopic-vorticity value a, 2) non-uniform 
steady states or coherent structures, 3) negative-temperature states (already 
predicted by Onsager's work on point vortices \T9\ ). Theoretical predictions 
were compared with numerical simulations and with experimental measure- 
ments in quasi-two dimensional fluids, e.g., in [T^l [HI 131 123 [!]■ However, 
under standard laboratory conditions fluids are viscous. Similarly, numeri- 
cal simulations require the introduction of a non-vanishing (hyper)viscosity 
in order to avoid some numerical instabilities and other artifacts. A non- 
vanishing viscosity can lead to noticeable effects, like the breakdown of the 
conservation laws, even at high Reynolds' numbers, i.e., when the viscosity 
"is small" . This is true especially when studying long-time processes like the 
generation of the coherent structures. In spite of the dissipative nature of 
the flows studied, in many cases, it was found that the agreement between 
the theoretical predictions based on the Miller, Robert and Sommeria (MRS) 
inviscid theory and (numerical) experiments was better than expected. 

The main purpose of the present paper is to better understand these is- 
sues in a more dynamical setting. We start by reviewing the MRS theory 
in Section [2 In Section [31 we consider viscous flows and propose a fam- 
ily of model evolution equations for the vorticity distribution (f){a,x,y,t). In 
Sections [IHSl we discuss a class of time- dependent distributions that maxi- 
mize a mixing index under certain constraints. In particular, in Section [S] 
it is shown that the time-dependent Lagrange multipliers appearing in these 
distributions can be chosen in such a way that the masses associated with 
each microscopic-vorticity value a are conserved. When these masses are 
conserved, the mixing index is the time-dependent analogue of the entropy 
functional used by MRS and it satisfies an H-theorem, as it is shown in Ap- 
pendix A. The distribution ^^(cr, ip{x, y)) associated with a given QSS can be 
obtained, at least in principle, by addressing the reconstruction problem, i.e., 
the question of how to extract its defining parameters from the QSS's u — ip 
relation. This is discussed in Section [3 In doing so we provide a natural 
framework for a time-dependent statistical theory connecting an appropriate 
initial distribution to the QSS distribution associated with the experimental 
uj — ip relation and evolving in agreement with the Navier-Stokes equation. 
The relation of the present work with the MRS theory and with the yard- 
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sticks' conditions of reference [0] is discussed in Section [HI 



2 The Miller-Robert-Sommeria (MRS) The- 
ory 

2.1 Review 

The pillars on which the statistical mechanical approach [TTlElj stands are the 
conserved quantities of the non-dissipative Euler equations. These quantities 
are: 1) the total energy Eg := fj^dxdyvl where 'v'o{x,y) is the initial 
velocity field and A is the area occupied by the fluid and 2) the area density, 
denoted by g{a) da, occupied by fluid with vorticity values between a and 
a + da, which is 



g[a) := / dxdy(f){a,x,y), 
J A 



where 0((T, x,y) > is the probability density of finding a microscopic vor- 
ticity value a at position {x,y), with 

da (j){a,x,y) = 1, for all {x,y). (1) 

For the sake of simplicity, we will ignore other conserved quantities which 
may be present due to some symmetries of the domain A, like the linear 
momentum and/or the angular momentum. Averages taken over this distri- 
bution will be indicated by (■■■), i.e.. 



(/(c^, X, y)) ■■= J da f{a, x, y)(f){a, x, y). (2) 

The macroscopic vorticity uj{x,y) is then 

u;{x,y) = (a) . 

It is assumed that the initial area density, denoted by goicr), is given by 

goia):= dxdyS{a-uJo{x,y)), (3) 



A 



where uJo{x,y) = V x is the initial macroscopic vorticity field. 
As usual, one derives the probability distribution for observing, on a mi- 
croscopic level, a vorticity value a by maximizing the entropy S under the 
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constraints defined by the conserved quantities. The entropy used by Miller, 
Robert and Sommeria (MRS) is 



S := —A ^ j dxdy J da (f){a,x,y)\n(l){a,x,y), 



(4) 



The probability distribution 05(cr, x,?/) that one obtains by maximizing the 
entropy S is 

(J^sicr, X, y) ■■= exp \-i5a^{x, y) + /i(cr)] , (5) 
with Z{'ijj{x, y)) := / da exp [—l3a'ijj{x, y) + /i(cr)] , 



where ^^(x, y) is the stream function, i.e., — A'?/'(x, y) = ujs{x, y), with ujs{x, y) 
the macroscopic vorticity field in the most probable state, i.e.. 



ujs{.x,y) = J daa(l)sia,x,y). 
We have then the u-^' relation 

ujs{x,y) = n{ilj{x,y)), 

with Q{ip) := J da a(f)s{cr,i'), 

where ^^(a, ip{x, y)) := ^^(a, x, y), 

which, in an experimental context, is often called the scatter-plot. This is a 
mean- field approximation, valid for this system [17] , see also eq. below. 
In eq. ([5]), /3 and fi{a) are Lagrange multipliers such that the QSS energy Es 

^s = \ I dxdyujs{x,y)ij{x,y) 



A 



2 

and the QSS microscopic-vorticity area distribution gs{<y) 



gs{a):= dxdy(f)s{,cr,x,y), (6) 



A 

have the same values as in the initial vorticity field, i.e., Es = E^ and 
gs{(y) = go{(^)- The system of equations is closed by 

^s{x,y) = -Aifj, (7) 
i.e., niip) = -Aij. (8) 
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which embodies the mean-field approximation. 

Since by construction gsio') = go{,o'), the assumption expressed by Eq. 
imphes that M„ the microscopic vorticity moments in the quasi-stationary 
state, 



M„:= jdaa^gs{a), (9) 
equal the macroscopic vorticity moments in the initial state, 

j daa'^goicr) = [ dx dy uj"^ {x , y) , 

confer equation ([3]). 



A 



2.2 The initial state 

As we have seen, in its standard formulation, the MRS theory introduces an 
asymmetry in the microscopic characterization of the initial state and that 
of the final quasi- stationary state. Namely, it is assumed that the initial 
state is an "unmixed state" for which the microscopic and macroscopic de- 
scription coincide. This allows to determine (?o(o") the microscopic vorticity 
distribution of the initial state from the macroscopic vorticity field uJo{x,y), 
confer eq. By contrast, the microscopic vorticity- area density of the most 
probable state, gs{o'), is given by fdxdy(f)s{cr,x,y), confer eq. In most 
cases, Gs{(j), the macroscopic vorticity density of the most probable state, 

Gs{cr) := / dxdy6{a - ujs{x,y)), (10) 
J A 

is different from the microscopic vorticity-area density of the most probable 
state, i.e., Gs{cr) 7^ gsio')- For example, for non- vanishing viscosity the even 
moments of 6*5(0") denoted by T2n, i-e., 

T2n-= / dxdyuj'^^{x,y), 

J A 

are smaller than those of gsio') since in the MRS approach one assumes that 
gsio') = go{o'), confer ( l25l) below. This difference between the moments of 
the macroscopic and those of the microscopic vorticity has led to confusion 
and to discussions in the literature related to the interpretation of the MRS 
theory [HIE]. 

According to the MRS theory the microscopic-vorticity moments M„ of the 
QSS, confer Equation ([9]), equal the microscopic-vorticity moments of the 
5-type initial distribution defined in (JSj). In Subsection III B of an earlier 
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paper [U] we have expressed this infinite set of equahties in terms of the uj — ip 
relation and the initial macroscopic vorticity field uJo{x, y). In order to assess 
the validity of the MRS approach, we introduced then a yardstick associated 
with each moment M„, see further the discussion in Section [9l 



3 Microscopic Viscous Models 

From here onwards we deal with time dependent probability densities, i.e., 
(f){a,x,y,t)da is the probability of finding at time t a microscopic vorticity 
value in the range (cr, a + da) at a position (x, y) which is non-negative and 
normalized 

(i(J 0((T, X, t) = 1. (11) 



The macroscopic vorticity field is 

u{x,y,t) = (a) , (12) 

where the notation introduced in ([2]) has been extended to the time-dependent 
case. In the inviscid case, the time evolution of 0((j, x, y, t) can be taken to 
be 

d(j){a,x,y,t) 

+ v[x,y,t) ■V(l)[(T,x,y,t) = 0, (13) 

where the macroscopic, incompressible velocity field v{x, y, t) satisfies ap- 
propriate boundary conditions and is related to the macroscopic vorticity 
^{x,y,t) hy 

V XV = ujz, (14) 

with z a unit vector perpendicular to the {x, i/)-plane. Consequently, the 
advective term in equation (fT3ll is quadratic in cf) and different values of a 
are coupled by this term. 

In this Section we introduce some model evolution equations for the prob- 
ability density (f){(7,x,y,t). We demand that these microscopic models be 
compatible with the macroscopic Navier-Stokes equations. The general form 
of these models is 

d(f>{a,x^, y, t) _^ _ V0(a, x, y, t) = vO{a, x, y, t), (15) 

with V the fluid viscosity and O as yet undeflned but constrained by 1) the 
conservation of the total probabilty J da (j){a, x, y, t), therefore. 



/ 



da = 0, (16) 
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and by 2) the macroscopic Navier-Stokes equation should follow from the 
microscopic model, therefore, 



j daaO = Auj{x,y,t), (17) 

so that, multiplying both sides of equation ( |T5l) by a, integrating them over 
a and making use of (|T^ and (fT7|) . one gets the Navier-Stokes equation 

duj{x^y,t) ^ _ \/u;{x, y, t) = uAu}{x, y, t). (18) 

For future convenience, let us introduce 

O =: A(f){a,x,y,t) + 00, 

with 

(O) = 0, (19a) 
and (aO) = 0, 

so that the constraints ( fT6|l and ( ITTI) are satisfied. 

It is convenient to introduce the "masses" m(cr, t) associated with each 
value a of the microscopic vorticity, 

m{a,t) := J dx dy (f){a, x,y,t). (20) 

In the inviscid case, incompressibility and the fact that the vorticity is just 
advected by the velocity field, allow for a complete association between a 
vorticity value a and the area that is occupied by such value. In this case 
one talks of 'conservation of the area occupied by a vorticity value'. In 
fact, in the inviscid case, z/ = 0, equation (|T5l) has a solution (j){a,x,y,t) = 
5(cr — uj{x, y, t)) so that the above-defined m(cr, t) is indeed the area occupied 
by the vorticity field with value cr. As soon as we introduce a diffusion process, 
as it is implied by equation ffT^ with u ^ 0, such an identification becomes 
problematic if not impossible. It is for this reason that we call m{a, t) a 
"mass" and, by doing so, we stress the obvious analogy with an advection- 
diffusion process of an infinite number of "chemical species" , one species for 
each value a. 

The time derivative of these masses is 



dm{a, t) 
di 



dxdyO(l){a,x,y,t). (21) 
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In order to derive these time-evolution equations fl^T]) it was assumed that 
there is no leakage of x, t) through the boundary, i.e., that the total 
flux of 0((J, X, t) through the boundary vanishes. 



dsn ■ 



v{x, y, t)(f){(T, X, y, t) - uV(t>{a, x, y, t) 



0, 



where the path of the integral is taken over the boundary and n is the normal 
unit vector. These conditions, as well as other boundary conditions that will 
be used in the sequel, are satisfied, e.g., in the case of periodic boundary 
conditions as well as by probability fields whose support stays away from the 
boundary at all times. In the case of impenetrable boundary conditions one 
has that, on the boundary, n ■ t? = 0, so that the last condition reduces to 

y rfs n ■ V0(a, X, y, t) = 0. (22) 

On a macroscopic level this implies, e.g., that 

ds n ■ Viu{x, y,t) = 



and the conservation of the total circulation J dx dy u!{x,y,t). 

The simplest model satisfying the above requirements is the one with 
= 0, i.e., 

d4>{a,x^, y, t) ^ _ V0(o-, x, y, t) = uA(j){a, x, y, t), (23) 

with the incompressible velocity field related to the vorticity as in eq. ( |T4l) . 
In spite of its simplicity, this model is very instructive because, while it 
dissipates energy, it has an infinite number of conserved quantities. Indeed, 
the masses are conserved, 

= 0, (24) 
dt ^ ^ 

confer (EI]). 

One of the consequences of the conservation laws is that all the 
microscopic-vorticity moments Mn{t) := Jdcr a^m{a,t) are constants of the 
motion, i.e., 

dMn 
dt 

In the sequel we shall assume that the conservation of all the microscopic 
moments M„ implies in turn that the masses m(cr, t) are conserved. This is 
the case if certain technical conditions are satisfied, see e.g., |26j . 
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On the other hand, all even moments of the macroscopic vorticity, call 
them T2nit), i.e., all the quantities 



^2nit) ■■= J dxdyu "(x,?/,t), 
in particular the macroscopic enstrophy J dxdyuj'^{x,y,t), are dissipated, 

dT2n 



dt 



-v2n{2n - I) j dx dy cu^^''-^^ \\/ujf < 0, (25) 



as it is implied by the Navier-Stokes equation ffTBl) . Also the energy E 
1/2 J dxdyv'^ is dissipated 



= I dxdyv-Av = —u / dxdy 



^v+r^v+2r^v 



dx'^ ) ^ \dy'^ ) ~^ \dxdy J 



<0, 



where 'ip{x,y,t) is the stream-function associated with v{x,y,t) and appro- 
priate boundary conditions, e.g., periodic boundary conditions, have been 
assumed. Under these boundary conditions, one also has, 

1 f 

= - dx dyuip, 
A dE f,,2 

and — — = —u dxdyuo . 
dt J 

Notice that the microscopic energy and the macroscopic one coincide, 

— - J dx dy J da J dx' dy' J da' a(l)(a, x,y,t)K{x,y, x' ,y')a'(l){a' , x' ,y' ,t) = 



- I dxdy I dx' dy' uj{x,y,t)K{x,y,x',y')uj{x',y',t) 



2 



IX ay iuyx,y, i)i\yx,y, x' 

'A J A 

In this last expression, K{x, y, x', y') is the Green function solving 



AK{x, y, x', y') = 6{x - x')6{y - y'), 

and the appropriate boundary conditions. 

In conclusion: The viscous Navier-Stokes equation l[T8\) does not exclude 
the possibility of an infinite number of conserved quantities m{a, t) as defined 
in equation (d^). 

In the following Sections we will consider more generic models with O ^ 
that conserve the masses m{a,t). 
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4 Chaotic mixing 



Considering once more the models defined by equations (l23ll and ( |T5l) . we 
notice that they are advection-diffusion equations for a non-passive scalar 
with the viscosity v playing the role of a diffusion coefficient. This type 
of equations has been studied extensively, see e.g. [20] and the references 
therein. It is well-known that a time-dependent velocity field {7(x, t) usually 
leads to chaotic trajectories, i.e., to the explosive growth of small-scale 0- 
gradients. These small-scale gradients are then rapidly smoothed out by 
diffusion, the net result being a very large effective diffusion coefficient, large 
in comparison to the molecular coefficient v. 

Based on these observations we will consider situations such that dur- 
ing a period of time, starting at to and ending at Ts when a quasi-stationary 
structure (QSS) is formed, a fast mixing process takes place leading to a prob- 
ability density, denoted by 05 (a, x, y) , that maximizes the spatial spreading 
or mixing of the masses m((j, Ts) and is such that ^^(a;, y) = Jda a(j)s{cr, x, y) 
satisfies the experimentally found uo — ip relation. More specifically, we will 
investigate solutions of equation ( |T5l) for suitable vO on the right-hand side 
satisfying the following conditions: i) at time t = Ts the solution 0((T, x, t) 
is the above-mentioned maximally mixed (f)s{o',x,y) complying with the ex- 
perimental uj — ip relation and ii) the mixing takes place much faster than 
the changes in the masses m{a,t), so that m{a,Ts) = m{a,tQ). In order to 
express these ideas in a quantitative form, we need a mathematical defini- 
tion of the degree of spreading or mixing of a solute's mass in a domain. As 
described in Appendix A, the 0-order degree of mixing so(cr, t) defined as 



SQ^a.t) := —A ^ dx dy (f){a,x,y,t)\n[A(f){a,x,y,t)/m{a,t)] , (26) 



is the right quantity in order to quantify the mixing of the microscopic- 
vorticity masses m(cr, t). Accordingly, the total vorticity mixing at time t, is 
measured by 



So{t) = —A ^ da dx dy (f){a, x,y,t) In (l){a, x,y,t)+ A ^ dam{a,t)\n[A ^m{a,t)^ . 



The first term coincides with the MRS entropy (jlj). If the masses are con- 
served, i.e., if m{a,t) = m (ex, to) , then the second term is constant in time. 
This allows us to give a mathematical definition of fast mixing, namely: fast 





More explicitly. 




(27) 
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mixing of m{a,t) takes place in a time interval {tQ,Ts) if the following in- 
equalities hold, 



dso{a,t) 



dt 



> 



A 



dm{a, t) 



dt 



for to<t< Ts. 



{2i 



Accordingly, the corresponding microscopic vorticity probability density 05(o", x, 
maximizes the total degree of mixing So(Ts) = J da so{a,Ts) under the con- 
straint that the change in the masses m{a, t) during the time interval [to, Ts] 
is very small. 

It turns out that, for our purposes, a second constraint is needed. Here 
we discuss two possible choices of this second constraint. One possible choice 
of the second constraint consists in using the macroscopic vorticity u^x, y, t), 
i.e., the solution of the Navier-Stokes equation with initial condition Uo{x, y) 
as input. This means that one demands. 



da 



ac 



{a,x,y,t) = u{x,y,t). 



Two points should be stressed: 1) since Lj{x,y,t) is a solution of the Navier- 
Stokes equation ffTSj) this choice of the second constraint is totally compatible 
with the presence of viscous dissipation and 2) this choice of the second 
constraint can be imposed at all times t, not only at time Ts when the QSS 
is present. The probability distributions which are obtained by maximizing 
the total degree of mixing 5*0 (t) under the constraints of given values m{a, t) 
for the masses (l20l) and given the vorticity field uj{x,y,t), are presented in 
the next Section. 

The second possible choice of the constraint stems from physical evidence 
that, in high Reynolds' number, two-dimensional flows, the energy is trans- 
ported from the small to the large scales and, consequently, it is weakly 
affected by viscous dissipation. More formally. 



dSoit) 



1 
dt 



dE 



dt 



for to<t<Ts 



(29) 



Therefore, one maximazes 5*0 under the constraints that the masses and the 
energy at time Ts have some given values m{a,Ts) and E{Ts). Introducing 
the corresponding Lagrange multipliers f3 and Jl{a), as well as a Lagrange 
multiplier 'y{x,y) associated with the normalization constraint ( ITTi) . the con- 
strained variation of the degree of mixing 5*0 leads to 



= paipix, y) + + 7(x, y) + In 



A(t)s{a,x,y) 



m[a] 



(30) 



I.e., 



0s(cr, x,y) = A m{a) exp{-l3a^p{x, y) - fi{a) - -/{x, y)), 
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where m(cr) := m{a,Ts). Since A^^m{a) > 0, we can define fi{a) := — /I(cr) + 
In A~^m{a) and implementing the normahzation constraint ([T]), one arrives 
at the probabihty density ([5]) of the inviscid, statistical mechanics approach. 
The new elements here are the conditions expressed by ( l28l) and ( l29l) . i.e., 
a criterium for the applicability of these equations in the case of viscous 
fiows. Notice, moreover, that in contraposition to the statistical mechanical 
approach, we do not require that the energy E(Ts) and the masses m{a,Ts) 
at time Ts be equal to their initial values; equations (j28ll and (!29l) only ex- 
press that the changes in these quantities are much smaller than the change 
in total vorticity mixing in the time interval {tQ,Ts). In the presence of vis- 
cosity E{Ts) < E(to), when the equality holds one recovers exactly the MRS 
expressions. 

In closing this Section, it is worthwhile recalling that under the physical 
conditions leading to the inequality in (l29l) the changes in energy are usually 
much smaller than the changes in enstrophy. This fact is at the basis of the 
so-called selective-decay hypothesis [IS], i.e., the conjecture that the QSSs 
correspond to macroscopic vorticity fields with energy E{Ts) which minimize 
the macroscopic enstrophy r2(r5). See also [Ij. 



5 The time-dependent extremal distributions 

In this Section we investigate the time-dependent probability distributions 
which are obtained by using the macroscopic vorticity uj{x,y,t) as second 
constraint. As it will be shown in the following Subsection, the form of these 
distributions is, 



The functions xi^, V, t) and /i(cT, t) will be called the "potentials" . In contrast 
to the situation in the statistical mechanics approach, these potentials can 
be time-dependent. Evidently, for these distributions one has 4>*{cr, x, y, t) =: 
(f){a, x{x, y, t),t), therefore the (x, y)— dependence of uj{x, y, t) as well as that 
of all the higher-order local moments mn{x, y,t) := (cr") and the centered lo- 
cal moments Kn{x, y, t) := {{a — lj)"') is only through x(a;, y, t), i.e., Kn{x, y, t) =: 



(p* (cr, x,y,t) = Z ^ exp [fi{a, t) + x{x, y, t)a] , 



(31) 



with 




(32) 
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Kn{x{x,y,t),t), moreover, 



uj{x,y,t) =: n{x{x,y,t),t), 

J da cr" exp [/u(cr, t) + xcr] 



and m„(x,t) := 



J da exp [fi{a, t) + xo"] 



moreover — = A2(x,^), and K3(x,t):= ^ • 

As it is easily checked, the local moments rhniXi^) satisfy then dfhn/dx = 
fhn+i — ^frin, i.e., the following recursion relation holds, 

fnn+iix,t) = Comnix,t), (33) 
d 

with C:= — + n{x,t). (34) 
dx 

5.1 Maximum mixing 

Time- dependent distributions as in Equation (IHTj) are obtained by maximiz- 
ing 5*0 (t), the total degree of mixing at time t, under the constraints of i) 
normalization, confer Equation ( iTTl) . ii) given values m{a^t) for the masses 
at time t as defined in ( I20l) and iii) given vorticity field Lj{x,y,t) at time t, 
i.e., the distribution's first moment (fT^ . (a) = uj(x,y,t). Indeed, introduc- 
ing the corresponding time-dependent Lagrange multipliers 7(x, y, t), Jl{a, t) 
and x{x,y,t) and denoting the extremal distribution by (j)M{o',x,y,t), one 
has that the vanishing of the first-order variation reads 

n { +^ I ^ +^ I +^ I 1 MM{(T,x,y,t) 

= 7(x, y, t) + x{x, y, t)a + /x(a, t) + In — , 

m((T, t) 

from which, in analogy with the derivation of the statistical mechanical for- 
mulas in Section m one obtains that 0A/(cr, x, y, t) has precisely the form 
of (f)*{a,x,y,t) in (EI]), i.e., (j)M{<y,x,y,t) = (j)*{a,x,y,t) with /i(cr, t) := 
— /I((T, t) -|-lnA~^m((T, t). The potential functions x{x,y,t) and ^{a,t) should 
be determined from the two following constraints 

^{x, y, t) = Z'^ jdaa exp [axix, y, t) + n{a, t)] , (35) 

and m(a,t) = J dx dy exp[ax{x,y,t) + fj.(a,t)] , 

where the vorticity field uj{x,y,t) is not a stationary solution of the Eu- 
ler equations but a time- dependent solution of the Navier-Stokes equations. 
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The main differences with the maximizer 05((T, a;,?/) of the previous Section 
are that these distributions (l)M{<^,x,y,t) are time-dependent and that now 
there is no energy constraint, therefore, the Lagrange multiplier x{x,y,t) is 
not a linear function of the stream-function ip{x,y^t). Up till now, the time 
dependence of /i(cr, t) has been arbitrary, in the next Subsection this time 
dependence will be determined such that the masses m((T, t) are conserved. 



5.2 The time-dependence of /i(cr, and the conserva- 
tion of the total moments 

Assume that at all times the probability density has the form given in ( |3T1) 
which is a time dependent probability density such that 5*0 (t) attains its 
maximum value compatible with (a) = uj{x, y, t) and given m{a, t). Inserting 
these expressions in the Navier-Stokes equation (fTSjl and making use of simple 
algebraic equalities like. 



d(f){a,x,y,t) 
dt 



(36) 



and 

VK2 = Ksix,t)Vx, 

with Kn '■= (((T — u;)") , leads to the following relationship between dx/dt 
and dfi/dt, 



dx 
dt 



a - u 



dt 



(37) 



Using this expression in order to eliminate dx/ dt from fl36l) . one finally arrives 
at Equation ffT^ with 0{a,x,y,t) expressed in terms of xi^^U^i) ^-^id of 
dn/dt. Namely, one obtains that 



0{a,x,y,t) = V ^ 
+ 



d[i 
'dt ~ 



dfi 
'dt 



V 



_i(fT - a;) 



d[i 
'dt 



+ 



K2 + ^ {a - uj) - {a - iuY 

J^2 



(38) 
(39) 



14 



As one can check, this expression satisfies (O) = and (cO) = indepen- 
dently of the specific form of dfi/dt. In other words, the (x, 2/)-dependent con- 
straints ( fTTl) and ( fT2|) hold at all times and for all possible time-dependences 

of /i((T, t). 

From ( l38i) it follows that the simplest viscous model of the type given 
in Equation (IHTI) . with an identically vanishing 0{a,x,y,t), can only be 
realized under very special, and rather trivial, conditions. Indeed, since 
0{a,x,y,t) = has to hold for any value of a, equation ( !38l) implies that 
djji/dt must be quadratic in a and that |Vx|^ may be time-dependent but 
must be (x, y)-independent. Therefore, the simplest viscous model with an 
extremal distribution as given in fl3T|) . 0{a,x,y,t) = and satisfying equa- 
tion flHHl) can hold only if x{x, y, t), and lj{x, y, t) are linear functions of the 
space coordinates. 

For general /i(o", t) there is no conservation of the masses m{a, t). However, 
choosing a suitable time- dependence of //(a, t) such that 



dxdy(j)0 = 0, 



(40) 



ensures the conservation of the masses, i.e., m{a,t) = m{a,to) at all times t, 
confer equation (12T1) . Using equation (l38il the last equality is seen to imply 
that, 



dfi((T,t) / dfi((T,t)\ (cr — Lu) I dfi(a,t) 

' (a-uj) 



dt 



dt 



dt 



(cr, x, t) dxdy = 
(41) 



K2 + — {a-uj)-[a -uj) 



|Vx|^0(o-,a;,|/,t) dxdy 



This equation is a complicated integro-differential equation for the time- 
dependence of /i(cr, t), however, using a Taylor expansion /i(cr, t) = fikit)cr'^, 
we can derive an infinite set of linear differential equations for the d^k/dt. 
In fact, multiply equation ( l38l) by cr"0(cr, x, ?/, t), integrate it over a and get 
then that: 



u {a 



-V 



fc=2 



dt 



with 



hnk{x, y, t) ■= ruk+n - mkirin - 



where m„ are the local moments mn{x, y,t) := (a") . The conservation of the 
moments M„ = J dxdy mn{x , y , t) requires then that Jdxdy Jdaa^cpO = 
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and hence that, 




dxdyhnk = ^ dxdyhn2 |Vx| , 



n = 2,3,.... 



(42) 



From this infinite set of equations, hnear in dfi2/dt, dfi^/dt, the dfik/dt can 
be solved, in principle. We have then a viscous model with an infinite number 
of conservation laws. Such a viscous model becomes physically more relevant 
by making it compatible with a quasi-stationary distribution (psi^', 'ipix, y)) as 
given by Equation ([5]) with associated relation at time Ts, i.e., at time 
Ts we can associate with 05 (a, i{^{x, y)) a distribution function as in Equation 
( 13T1) which is a solution of the time-evolution Equations ( ITSl) and ( |38l) . with 
suitable initial conditions such that (j)*{a,x,y,Ts) = (p{a,x{x,y,Ts),Ts) = 
(ps{cr,ipix,y)) with x{x,y,Ts) = -Pip{x,y) and with^/x(cT, T5) = n{a). The 
question of reconstructing the distribution function (j)s{o',ip{x,y)) from the 
experimental u — ip relation will be addressed in Section [71 

In concluding we want to remark that, as it is shown in Appendix A, 
an H-theorem holds for the extremal distributions (f)*{cT, x, y, t) in the case of 
conserved masses m{a, t) treated in this paper. More precisely, it is for these 
distributions with conserved masses that one can prove that dSQ{t)/dt > 0. 
Such an H-theorem does not hold for the other measures of mixing Sj. with 
r 7^ which are defined and discussed in Appendix A. Therefore, this result 
fixes the 0-degree of mixing 5*0 defined in (l26l)-(l271) as the appropiate measure 
of mixing. 

6 Fast-mixing condition 

In this paper we are mainly concerned with dynamical models such that 
dm{a,t)/dt = 0, as treated in the previous Subsection. In such a case, the 
inequalities in fl^^ will always be satisfiecl^ and assuming that the probability 
distribution is of the extremal form given in Equation ( !3T|) . we can make use 
of equation ( 1671) in order to write inequality ( l29l) of the fast mixing condition 
as. 



Two observations are in place: the viscosity drops out from the final expres- 
sion and the r.h.s. is totally determined by the macroscopic vorticity field 
uj{x,y,t). 

^We exclude the exceptional situation dso/dt = 0. 




dxdy (|Vln0|^) ^ j dxdyuj"^ 




(43) 
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Although the energy-related fast-mixing condition fl29|) does not play a 
role in obtaining the time-dependent distributions 4>m{o', x, y, t) = (f)*{a, x, y, t), 
it is still an interesting issue to determine to which extent this fast-mixing 
condition is actually satisfied or not by the 0*(cr, x, y, t) distributions or more 
general 0(cr, x, t) satisfying equations (fT5l)-( |T7|) . In fact, for high Reynolds' 
number one expects that after a time-interval of fast mixing one arrives at 
a well-mixed probability distribution with E{Ts) < E{to). In order to in- 
vestigate this in more detail we will 1) derive a lower bound to (|Vln(/)|^) 
and 2) we will investigate the extrema of the l.h.s. of equation ( H3i) taking 
into account the constraints ( JTTl) and ( fT2l) . As it will be seen, the special 
distributions defined in equation fl3T|) play a mayor role in both questions. 



6.1 Lower bound on {|||Vln0|^ 

In order to find a lower bound for (|Vln0|^) in terms of |Vci;|^ we begin by 
noting that Vcu = {aVliKp) can also be written as Vcu = ((cr — cj) V In 0) 
because (Vln0) = 0. Applying then the Cauchy-Schwartz inequality to 
|Vti;|^ = |((cr — u;) V In 0)1^ leads to the desired lower bound, 

< {{cx-ujf){\V\n(j)f) =K2(|Vln0|^), (44) 
i.e., ^<(|Vln0r>, 

where K2 is the centered second moment K2{x, y, t) := j da (a — u!{x, y, t))^ (f){a, x, y, t). 

The lower bound on {||Vln0|^) that we have just found means that the 
fast-mixing condition fl43p holds whenever 

1 /■ , , \Vuj\^ 1 



^ J dx dy ' ^ 'e j ^^^y^"^' 

irrespectively of all the higher-order moments with n > 2. Recall that J dx dy | Vcul 
is directly related to the dissipation of the macroscopic enstrophy confer 
equation fl2^ with n = 1. 

It is rather straightforward to determine the family of vorticity distribu- 
tions, let us call them 4>lb{<7, x, y, t), that reach the lower bound in ( H4l) . i.e., 
that satisfy (|Vln02.B|^) = -^2""*^ l^'^l^ • Appendix [B] it is proved that 

0LB(a,X,?/,t) = 0*((T,X,?/,t), 

with (j)*{a, X, y, t) as given by equation (1311) . In the present context, the input 
for the determination of the potential functions fi{cr,t) and xi^iVi'^) ^-^e the 
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first and second a— moments, 

uj{x, y,t) = j da acpLsio; x, y, t), (45) 

and K2{x,y,t)= da {a - uj{x,y,t)f (f)LBicr, x,y,t). 



6.2 Extrema of Jdxdy ( |Vln0| 



In Appendix O, we investigate the extrema of J dxdy (|Vln0|^) taking into 
account the (x, |/)-dependent constraints (fT^ and (fTTl) . It turns out that in 
order to obtain sensible solutions it is necessary to constrain also the distri- 
bution's second moment (cr^) . We find that all the extremizer distributions, 
call them (j)ext{o',x,y,t), are local minima of Jdxdy (|Vln0|^) and have 
the form 4>ext{o', x, y, t) = (f)*{a, x, y, t), with (f)*{a, x, y, t) as in equation (I3T11 . 
This is in total agreement with the lower bound (^IV liKpiBf) found in the 
previous Subsection, confer Appendix [Bl In fact, as one can check, 

i.e., these extremal distributions (pext reach the lower bound (Hlj). In the 
present case, the potential functions /i(cr, t) and y, t) should be deter- 
mined from the input functions uj{x, y, t) and K2{x, y, t), just as in equations 

Summing up the results of the two last Subsections: 1) all distributions 
(f){a,x,y,t) with first moment uj{x,y,t) and second moment m2{x,y,t) = 
iu'^{x,y,t) + K2ix,y,t), satisfy (|Vln0|^) > (iVln^isO = |Vcj|^ , 

2) the lowest possible value of (|Vln0|^) is achieved for (f)LB{o',x,y,t) = 
(f)ext{o',x,y,t) = (f)*{a,x,y,t), with (j)*{a,x,y,t) as in equation ( |3T1) and 

3) if Jdxdy {\Vujf /K2) > (A/E) J dx dy , then the fast -mixing condition 
( 143|) holds implying that the energy changes in the time interval [to, Ts] are 
small. 



7 Reconstructing //(a) from experimental data 

Suppose that in an experiment one is given an initial vorticity field with 
its corresponding energy Eo and that one finds, from a time T5 onwards, a 
quasi-stationary vorticity field 

^s{x,y) = n{^{x,y)), 
with ijjmin < i^{x,y) < i/j max; 
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with a monotonic and an energy Es < E^. As we show below this 

experimental data can be used in order to determine the potential 
occuring in the quasi-stationary distribution ^^(cr, ^Z'). Once the distribu- 
tion (psicjip) has been reconstructed from the experimental data we can 
then associate with it a time-dependent distribution function (j)*{a,x,y,t) 
as given by equation flHTl) which is a solution of the time-evolution Equa- 
tion ffT^ with suitable initial conditions and such that at time t = Ts one 
has (f)*{(r,x,y,Ts) = (j)M{(^,xix,y,Ts),Ts) = (f)sicr,i/j{x,y)) with x = -/^^ 
and with fi{a,Ts) = fJ'io^)- Here /3 can be defined such that M2(Ts), the 
second microscopic-vorticity moment of the QSS at time Ts, is equal to the 
enstrophy of the initial vorticity field u}{x,y,to). More explicitly, the con- 
dition M2(Ts) = M2(to) = T2 leads, also in the viscous case, to 

_ J^dxdy {dn/djj) 
'-2 '-2 

see further the discussion in Section 

The time-dependence of fi{a,t) is chosen as in Subsection 135.21 i.e., such 
that all masses m{a,t) are constant in time. 

In order to determine the potential from the experimental ^{ip), we 
notice that any monotonic function, like ^{x), can be expressed as 

nix) = ^inz{x), 

dX 



/oo 
da exp + xcr] 
-00 



for some appropriate function fi{a). From the above formula and demanding 
that X = —Pi', it follows that it is possible to express Z{x) in terms of ^{4'), 

Z(x) = exp[-/5 r di^'niij')] 

where (3 is determined from the experimental ^{ip) relation using equation 
( H6l) . In order to illustrate this procedure, we now treat a number of cases in 
which we can obtain /i(cr) from VL{%1)) either analytically or approximately. 

7.1 Linear u — ijj relation 

The case of a linear uj — ip scatter plot, i.e., 

or ni{x) = -aiX/P, 
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can be treated rather straightforwardly. In this case one gets, 

Ziix) = exp 

with Xmin •= —Pi^min- As One can check, taking /U/(cr) such that 
. . 1 



Q^l ( 2 _ 2 \ 
2^ \^ '^min/ 



exp iJ,i[a) 



exp 



«1_ 2 

2^-^111111 



exp 



2ai 



with ai/P < 0, leads to the desired Zi{x), 

/oo 
da exp [/xz(a) + xcr] = exp 
-oo 





2 






exp 


2^-^min 


exp 


2/3^ 



Since 



< 0, and 

Es = dxdyus'ips 



1 

"1 



dx dyuj\ > 0, 



it follows that ai is positive and /3 is negative. The corresponding probability 
density is then 



(j)s{cr,x,y) 
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/3 



exp 



^ 'a - al^p{x,y)f 



2ai 



i.e., a Gaussian centered on aiijj with a width which has the form of 

(IHl) with fi{a) = /i20"^ and fi2 = P/ (2ai) . In the present case the expression 
( H6l) for /? reads, 

= ^ (47) 

where T2it) is the enstrophy J dxdyLj'^{x,y,t). 

7.2 General nonlinear uj — ip relation 

In the more general nonlinear case we first notice that 
djji^a) 



da ■ 



da 



exp [ii{a) + X0-] = -xZ{,x) + exp [ii{a) + xo"]! 



+ 00 

— oo 



The boundary terms vanish, so that, using the notation of equation (j2]), one 
has 

'''^\ = (48) 
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Introducing into this equality the Taylor expansion of fi{cr), 

k=2 k=2 

where the coefficients /ifc are independent of x, one gets, 

kfikmk-i{x) = -X- (49) 

k=2 

Using the recursion operator C defined by fl33l) . this equation can be rewritten 
as, 

'fc-2 



fi = -X, (50) 

.fc=2 

where := o £. Equation flHIl) is a nonlinear differential equation in 
fi(x) but it can be reduced to a linear equation in the partition function 
Z{x)- This simplification is achieved by multiplying both sides of equation 
fl50|) by the partition function Z{x), and noticing that 

forn>l, Z(x)/:"°^^= ^ (51) 

and for n = 0, ZM^ = —■ 

dx 



One realizes then that (15011 reduces to the following linear differential equa- 
tion, 

5Z^^^^3rr = -^^' (52) 

k=2 ^ 

that the partition function Z{x) must satisfy. In Subsection 186.31 and in 
Appendix [D] we consider some cases in which Equation (I52j) can be reduced 
to a finite-order differential equation. 

7.3 Slightly nonlinear u — i/j relation 

It is often experimentally found that the u — ip plots satisfy ~ —Q{ip). 

The odd character of implies that /i(— o") = /u((t), i.e., 

yu(cr) = /i20-^ + /i4(j^ + fi^a^ H . 

Moreover, in many cases these plots are nearly linear. More specifically, in 
the Taylor expansion for yu(o") and in 

^{x) = fiX + f3x' + hx' + -; for |xl<Xmax (53) 
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one has that in a substantial interval around = or x = 0, 
Zix) = Cexp 



(54) 



with 



for n = 1, 3, 5, 



\fn+2\X^ < \fn 

and /i„+2 < /i2yUn for = 2,4, 6 



Inserting these powers expansions of Z{x) and fi{cr) into ( l52l) allows us to 
express the {/in} in terms of the {/«} , i.e., it allows us to determine the prob- 
ability density exp ^{a) from the experimentally known scatter-plot ^{x)- 
For example, retaining terms up to /s and /| in the Taylor expansion of 



Equation 



one gets that 

/i4 



— -^/s/i 

D 

/^8 = A'lO = 



j/3/1 
= 0. 
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8 Conservation of a finite number of total 
moments 

As we have seen in Subsection 135.21 the masses m{a, t) are conserved if the 
time dependence of t) satisfies equation (HTj) . For the sake of simplic- 
ity one can demand less, for example a finite number of total moments M„ 
for 2 < n < will be conserved by requiring that fi{<J,t) = J2iL2 f^ii^)'^^ + 
J^i^N+i /^i^* where the /ij with i > N+1 are time independent and dfi2/dt, . . . , dfi^/dt 
satisfy equation (H2l) . In this case the derivation of an H-theorem as given 
in Appendix |X] does not go through even if one assumes that /ij = 0, for 
z > A^ + 1, and consequently that J dadxdy (f)0 Incp A = holds, confer flU7|) . 
This is so because the m(cr, t) are not conserved and instead of equation flUTl) 
one has that 

^ = ^ / i I + i / da^-^ (ln,„(.) + 1) . 

However, considering that fast mixing takes place at times t <Ts for which 
the changes in masses can be assumed to be small, we still may use equation 
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as the fast-mixing condition. In the present case inequahty 
satisfied automatically but only in the following weaker version 



is not 



d_ 
dt 



da (T'^So{a,t) 



> 



A 



dMn{t) 



dt 



for n>N + l. 



8.1 Gaussian distributions conserving the second mo- 
ment 

It is instructive to consider the N = 2 case, i.e., 



dt 



dt 



-a 



Then from (1421) one gets that 

d^^2{t) _ Jdxdy [K^'Kl + Kl - K,] \Vx\' 



dt 



V- 



J dxdy [K^^Kl + Kl - 



If further fik = for k > 2 and /X2 < 0, we have then that In 
the corresponding solution is a Gaussian distribution. 



(55) 
and 



0G(o-,a;,2/,t) 



l/i2(t)| 



TT 



exp 



X{x,y,t) 

2/i2(t) 



In this case we have a linear uj — ip relation and since (a) = uj{x,y,t), it 
follows that x(x, y, t) = —2n2{t)uj{x, y, t), and K2{x, y,t) = — [2fi2{t)]^^ .The 
centered moments Kn{x, y, t) of these Gaussian distributions are independent 
of {x,y), therefore the differential equation fl3^ satisfied by fi2it) reduces to. 



dt 



~A j ^^^^ '^^'^ 
lJi\{t) j dxdy \Vuj\ 



(56) 



4z/ 



The solution is. 



^^2{t) 



/i2(0) 



[1 - 2A-i/i2(0)A(t)] 



(57) 



where A(t) := 2v j ds j dxdy \Vuj\^ > 0. 



From equation ( |25l) with n = 1, one sees that A {t) is the accumulated change 
in macroscopic enstrophy r2(t), i.e., A(t) = (r2(0) — r2(t)) . Since A (t) > 
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grows with time and /i2(0) < 0, it follows that |/U2(^)| decreases, i.e., the 
width 1/ \/JJi2(J)\ gets larger as time increases. 

If the initial distribution is sufficiently narrow so that |/i2(0)| ^ A [2A(oo)]^^ 
then equation (|57j) tells us that there must exist a time T' such that for all 
later times t > T' 



for t > r, /i2(t) ~ -A [2A(t)] 



-1 



i.e., the inverse width fi2(t) becomes independent of its initial value /i2(0). 
Often A(cx)) ~ r2(0) and the width will become independent of its initial 
value if |/i2(0)| ^ A [2r2(0)]^^ , which depends only on the initial conditions. 
In particular, for the 5— type initial conditions as in ([3]) one has that 



lim /i2(t) 

tio 



and Equation (1571) becomes 



The condition for fast mixing in the present case reads. 



[A-2/X2(0)A(t)]7 ' E 

If the initial distribution is sufficiently narrow, i.e., |/i2(0)| ^ [A(cx3)]~^ then 
the condition for fast mixing at t > T' becomes, 

2z/ /• , , ,2 A(t) /• , , 2 
— / dxdy \Vuj\ ^ ^, , / dxdyu . 
A J ^ ^ ' E{t) J 

2 

If the initial flow is unstable then |Vci;| can become very large; in fact, 
these gradients will become larger the smaller the viscosity u. Accordingly 
one conjuctures that even in the limit of vanishingly small viscosity one may 
have that 



2u 



lim — / dxdy I VcjI^ > 0, 



see the analysis presented in [lOj . 

The corresponding dynamical picture would be as follows: due to chaotic 
mixing A{t) grows very quickly up to a time T5 when uj{x, y, T5) becomes, to 
a good approximation, a function of the corresponding stream-function, i.e., 
u{x, y, Ts) — fl{ip{x, y, T5)). At this moment one has a stationary solution of 
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the inviscid Euler equations, the (weak) time dependence is due to viscous ef- 
fects. If the uj — ip relation is hnear then the quasi- stationary state (ps (o", a;, y) 
is in total agreement with the Gaussian solutions 0g (c", y, Ts) provided 
that the initial distribution is sharp enough, i.e., that |/^2(0)| ^ A [2A(oo)]~^ 
and that Ts > T' . In fact, from the linear u — ip relation uj = aiip of Section 
[7]A. with /3 satisfying fHUj) it follows that, 

= -r iftix = -A^' (r2(0) - T^iTs)) X = (2/i2(Ts))-' X- 

The Gaussian distributions that we have discussed above conserve only 
the distribtuion's second moment. Suppose that at a certain instant T we 
have a Gaussian distribution, i.e., that ^{a,T) = fj,20'^, can this distribution 
conserve all moments and stay Gaussian at later times t > Tl Equation ( HTl) 
tells us that for this to be true the following identity should hold. 



dt 



UJ 



(f){a, X, y, t) dxdy 



-V j [K2- {o - ujf]\Vx? 'P{(^.x,y,t)dxdy 



or, after some manipulations, that 
dn,2 



j [(cr - uof - K2] 0(cT, X, y, t) dxdy = 
j [(cr - ujf - K2] \Vxf 0(0-, X, y, t) dxdy = 
Therefore |Vx|^ must be independent of {x,y) and so must be |Vci;|^ since 



dt 

V 



X = 2fi2UJ. This corresponds to the simplest dynamical model, equation (123|) 
with = 0. Therefore, a Gaussian distribution can conserve the masses 
m((T, t) only in the very special case of a space independent | Vc(j|^ . Moreover, 
the last equation implies that, 

dfJ'2 1^ ,2 

= 4z//i^(t) \Vujf , 
in agreement with our previous results, confer Equation ( l56l) . 

8.2 More general distributions with conserved second 
moment 

Here we study the family of nonlinear /i(cr, t) with only one time-dependent 
(T-scale g(t), i.e., 

/i(a,t) = Jl{q{t)a), 



25 



and accordingly, 



dt 



(59) 



a dfi dq 
q da dt' 

where Jl' := djl{x) / dx\^^g^ and the time-dependence of q(t) will be deter- 
mined in the following paragraphs such that the second global moment is 
conserved, i.e., M2(t) = Jdxdy (a^) = M2(0). 

The recursion operator C defined in ( 1331) allows us to write 



dfj.\ 
'Taj 



\ da 



confer Equation ( HSl) . After some algebra, it follows from ( !59l) that, 

1 dq 




and i/((T^r) 



3 (1 + x^) 



dn 



3fi^ + x^"^ 



Multiplying equation flH^ by o"^, taking the average over the vorticity distri- 
bution and using the expressions we just derived one arrives at, 

1 dq 



y{a^O) = v\Vx[ 



K2 



q dt 



2Ko 



K2 



The conservation of the second global moment M2 = Jdxdy (a^) requires 
that f dxdy (o"^0) = 0, it follows then that M2{t) is conserved if, 

Idq _ Jdxdy \Vxf [Kj + Kl - 7^27^4] 
'qltt ~ ^ 7rf^rf^K2^^]2^|^^7l^^ ■ 

In most cases one will have that [Kl + Kl — K2K4] < while [2/^1 — ^IK^] > 
0, therefore, the inverse width will decrease in time, i.e., the width of the 
distribution will grow in time. If the distribution is Gaussian then K^ = 0, 
K4 = 3Kl , K 



— [2fi2{t)] ^ and the last expression becomes, 
ldq_ Jdxdy \VxfKl 



q dt 



J dx dy K2 



K2 

I 

A 



dxdy iVxT 
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This is in agreement with the result given in equation fl3Bl) . 

If djl{x)/dx is a nonhnear function of x then it does not follow that 
J dadxdy (f)0 In {(f) A) = and there is no H-theorem, confer Equation ( !66l) . 
Assume now that at time Ts there is a nonlinear u-^ relation which can 
be associated with a distribution of type ( 13T1) with a t) as in equation 
(13^ . We can choose the value of /5 according to (HUj) so that ^2(0) = 0. In the 
context of our non-Gaussian distributions with conserved second moment this 
means that the distribution (f){a,x,y,Ts) = 05(cr, ■?/'(a;, y)) can be obtained 
from a 5-function like initial distribution. 



8.3 Non-Gaussian distributions with conserved second 
and fourth moments 

Let us now consider the special case when in the Taylor expansion of 
only /i2 and /i4 do not vanish. Then the differential equation ( 152|) satisfied 
by the corresponding Z (x) reduces to, 

dZ d^Z 

2/X2--1 — h 4/X4--— J = -xz with < 0, (60) 
(^X "'X 

From this equation one determines the asymptotic behaviour of ^{x) = 
dlnZ/dx as 

^^%-»3 V4|fniy 61,^41141*141/ \ X J 

We see that in this case limf2(x) ±cxd for x ~^ ±00. In fact, from Equation 
fl3^ one can prove that ^{x) remains finite in the limit |x| 00 iff infinitely 
many do not vanish. Using equation ( I6OI) one can express the coefficients 
/i2 and /i4 in terms of the experimentally determined coefficients of ^{x) = 

fiX + /sX^ + fbX^ + Inserting this expression in equation fl60l) . expanding 

in powers of x and equating the coefficients of and x^ on both sides of the 
resulting equality leads to 

2/3 + 24/1/3 + 40/5 
4/4 + 36/2/3 + 8O/1/5- 24/1' 
h 



m + 36/2/3 + 8O/1/5 - 24/. 
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In the special case /X2 = the solution of equation (|60il with boundary 
conditions Z(0) = 1, Z'(0) = and Z"(0) = is the Meijer G-function or 



27 



hypergeometric function o-^2([| , f], |/i4|). The corresponding is 

plotted in Fig. 1. As expected for this special case with = 0, around 
X = the function fl{x) is not linear but cubic, i.e., /i = 0. 

The coefficients /X2 and /i4 in the case of equation (l60l) can be time- 
dependent; their time-dependence can be fixed, e.g., by requiring that the 
second and fourth global moments M2 and M4 be conserved. Such require- 
ment, i.e., equations (H^ with n = 2, 4 and /i„ = for n 7^ 2, 4 leads to 



~dt 



H22 J dxdy \Vxf h^2- H42 J dxdy {Vxf ^^22 

^22^44 ~ ^24^42 



dfX2 f dxdy |Vx| /122 d^iH2i 

= u- , 

dt H22 dt H22 ' 

where Hnm ■= / dxdy Km 



m=0 



m 



and the relation m„ = X] ( ) ^'"-^n-m so that the time- derivatives can 



be expressed as complicated functions of K2,K2,, ....,K8. It is instructive to 
consider a nearly Gaussian situation in which the Kn are replaced by their 
Gaussian values K2m~i = and K2m = (2m — Then 

h22 = 2Kl K2 = l2Kl{K2 + cj'), 

/i44 = 24 {Kl + lOfsTscu^ + 3cu^) , 
dfii _u f dxdy \Vujf (u;^ _ raA^^) i^a"^ 
"dT " 2 AKI + AK2T2 + 3{T4-Tl) ' 



y dxdy \\/ujf -6{K2 + TaA-^) 



AK^ J ^ ' ' V ^ ' ^ J ■ 

The expression for dii^/dt is governed by the correlation between |Va;|^ and 
uj'^. Close to the maxima of uj the gradient will be very small, analogously 
the gradient maybe large in the areas where uj is not close to an extremum. 
Therefore, it is reasonable to assume that this correlation will be negative, 
so that dni/dt < while dfj,2/dt > 0. This means that, in the context of the 
models with conserved second and fourth moments, a slightly non-Gaussian 
initial distribution will evolve into a more strongly non-Gaussian distribution. 



9 Discussion and conclusions 

In the MRS approach the conservation laws of 2D inviscid fiows play a cen- 
tral role. However, the unavoidable viscosity makes these conservation laws 
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invalid; in fact, at high Reynolds' numbers it is only the energy that is 
approximately conserved [15j. On the other hand, as we have seen, the con- 
servation of the microscopic- vorticity masses m(cr, t) is compatible with a 
non-vanishing viscosity, it leads to spatial mixing of the masses without de- 
stroying or creating them. This is the picture adopted in the present paper. 

By working with the family of maximally mixed states studied in Sec- 
tions miH] we reached a mathematical formulation that, as far as the QSS is 
concerned, shows many parallels with the MRS formulation. In Section [7] we 
addressed the problem of how to determine the QSS distribution 4>s{cr, x, y) 
from an experimental u — ip relation and j3 given by Equation ( H6l) . Identi- 
fying this (f)s{o',x,y) with the distribution (f){a,x,y,Ts) = (f)M{cr,x,y,Ts) = 
'pLB{<^,x,y,Ts) = (t>*{a,x,y,Ts) of Section and using a time-dependent 
fL{cr,t) satisfying Equations fHT]) - fH21) and such that fi{a,Ts) = ft^cr) we have 
a dynamical model that conserves the masses, i.e. with m{a,t) = m{a,to), 
and connects the experimental u — ip relation found at time Ts with an initial 
condition (j){a,x,y,to)- An extra bonus that follows from this methodology 
is that an H-theorem holds, moreover, it holds for one measure of spatial 
mixing, namely for the 5*0 degree of mixing, confer Appendix A. When the 
masses m{a, t) are conserved, this mixing measure 5*0 coincides with the MRS 
entropy. 

In Subsection III B of an earlier paper [6] we expressed the quantities 5„ 

5n := j da j dxdy [a" - ^^{x, y)] (psicr, x, y), 

in terms of spatial integrals of certain polynomials in and its derivatives 
d'^Q/dip^. For example, for n = 2, one has that 

62 = — — / dxdy ■ 



j3 Ja dip 

and the MRS assumption = ^2 + Ff — F^, can be used, as we did in Section 
[71 in order to determine (3 from the experimentally accessible quantities ^l{ip) 
and the initial and final enstrophies F2 and Ff as in Equation (H6ll . From 
these 6n one finds the corresponding initial values 

S° := jda j dxdy [cr" - y)] (f){a, x, y, to), 



because the conservation of the moments M„ implies that 5° = 5„ + F,^ — F 
where F,;^ := fdxdyu!g{x,y) and F° := f dx dy {x , y) are the macroscopic 
vorticity moments in the QSS and in the initial state, repectively. In the 
MRS approach the initial distribution is assumed to be a 5-function as in 
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Equation Q, consequently, all 5° should vanish. This led us in [6J to propose 
the quantities, 

n = 2,3,.... (62) 

n n 

as yardsticks in order to investigate the validity of the statistical mechanical 
approach according to which all these quantities should equal 1. Choosing (3 
as in P6l) the yardstick relation 82/ (F^ — Ff) = 1 is automatically satisfied 
while the remaining yardsticks for n > 2 are nontrivial checks for the valid- 
ity of the statistical mechanics approach. When all the yardsticks relations 
5nl (r° — r„) = 1 hold the quasi-stationary state predicted by the MRS ap- 
proach is in agreement with the experimental u — ip relation, moreover, it is 
also the solution of equation f|T5|) at time T5 with conserved total moments 
Mn and starting from the initial condition (l){a,x,y,to) = 8{a — uJo{x,y)). 
When not all the yardsticks relations are satisfied then the MRS approach 
gives an incorrect but approximate prediction of the experimental re- 
lation. In such it still may be possible to obtain a rather sharp initial 
distribution, which is determined by the experimental relation, such 
that the quasi-stationary state is the solution of equation ( IT^ at time T5 
and with conserved total moments Mn{Ts) = Mn{to)- 

Hence it maybe concluded that our class of dynamical models with non- 
vanishing viscosity and with conserved masses yields a self-consistent dynam- 
ical description containing the MRS approach with 5— type initial conditions 
in the case that the experimental values of the yardsticks relations are 1 and 
with different initial conditions in all other cases. In these other cases the 
predictive power of the statistical mechanics approach, i.e., the prediction of 
the correct QSS on the basis of a 5— type initial condition, is lost and the 
reconstruction of the correct initial distribution in the context of the models 
with conserved masses is a rather difficult if not forbidding problem. In such 
cases it may be more convenient to give up the conservation of all masses, 
require that only the second moment M2 be conserved and use the Gaussian 
distributions of Subsection I8B.1I if the u — ip relation is linear or the family 
of /i(cr, t) with only one time-dependent scale q(t) of Subsection 186.21 in the 
case of a nonlinear uj — ip relation. 
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A Degrees of mixing 



A mass of solute m achieves the highest degree of mixing when it is homoge- 
neously distributed over the area A, i.e., when its concentration is constant 
and equal to m/A. Let us introduce S{a, x, y, t) the spatial distribution of the 
cr-species at time t, i.e., 

r/ .^ <Pi(^,x,y,t) 

6{a,x,y,t) := > 0, 

m(cr, t) 

dxdy6{a,x,y,t) = 1. 

In order to determine how well mixed is this a-mass one has to determine 
'how close' is the corresponding 6{a,x,y,t) to the homogeneous distribution 
1/A. As it is knownp], given two spatial densities, called them Si{x,y) and 
62{x,y), there is a one-parameter family of non- negative, convex functionals 
dr{Si,S2) satisfying dr{Si,S2) = < — > 6i =' 62 and measuring a sort of 
distance between them, namely 

dri5i,52) := [r{l-r)]-^ jdxdy5i{x,y) [1 - {52{x,y)/5i{x,y)y] > 0, < r < 1, 
and, by imposing continuity in r, 

c^o('^i,^2) ■=- j dxdy6i{x,y)\n{62{x,y)/6i{x,y)) , 

di{6i,62) := - / dxdy62{x,y)\n{6i{x,y)/62{x,y)) . 



For r > 1 and for r < 0, they may diverge. However, since in our case 
62{x,y) = 1/A 7^ 0, values r < would also be possible. The dr are not 
distances because, e.g., for r 7^ 1/2, they are asymmetric in their arguments, 
i.e., dr{5i,52) ^ dr{52,5i). 

Accordingly, we can measure the mixing degree of the cr-mass by —dr{S{a, x, y, t), l/A). 
Its weighted contribution to the total mixing degree will be denoted by 
Sr{a,t), i.e., 

Sr{a,t) := m{a,t) [Ar{r - ^ dx dy 5{a, x,y,t) [1 - {1/A5{a, x,y,t)Y] < 0, for < r 
So{cr,t) := —m{(7,t)A~^ / dx dy 6{(7, x,y,t)\nA5{a, x^y,t) 



A(f){a,x,y,t) 



■m{a, t) 

Si((T, t) := m(cr, t)A~^ j dxdy hiA5{(T, x, y, t) < 0, 



<0, 
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the corresponding total r-order degree of mixing is then, 



Sr{t) := [Ar{r - 1)] ^ J da J dx dy (f){a, x, y, t) [1 - {m{a, t) /A(j){a, x, y, t)Y] < 0, 

(63) 

SQ{t) :— —A~^ J da J dx dy (p{a, x, y, t) In [^0(cr, x, y, t)/m{a, t)] < 0. 

Under certain conditions, these degrees of mixing satisfy a kind of H- 
theorem. In order to see this, let us compute the time derivative of Sr{a,t), 
with < r < 1, 



< r < 1, 



dsr{a, t) 
di 



Ar (r — 1) 



dx dy 



1 + 



(r — 1) m'^ 

{MY 



(90 r / m dm 
'dt~A\A^) 'dt 



and consider the mass-consevation models with dm/dt = so that dsr{a, t) / dt 
is totally determined by the first term in the curly brackets and one has that. 



dsr{a, t) 
di 



r (r — 1) A 

u f 
= - / dxdy- 



dx dy 



1 + 



(r — 1) m^{a) 



[-V ■ V0 + uA(f) + uO(f)] 
(64) 



(MY 



(f) r 



with < r < 1, 



While for r = 0, 



(65) 



dso{a,t) V 



dt 



^ Jdxdycf) 



V0 



dx dy I 





V0 



+ 01n.^'"' 



A(l) 



-01nA(/) 



(66) 



In the simplest viscous model with O — 0, the last formulas lead to a 
H-theorem for each cr- value. 



dsr(a,t) /■ , , m'(a) iVchl ^ , ^ 



In the more general viscous model with O ^ and conservation of the total 
moments, i.e., with Jdxdy (^a"'0') = 0, one can prove an H-theorem only 
for the total 0-th order mixing So{t). This is done by assuming that 4> is 



33 



of the form given in flHTl) so that ln0 can be expressed as a power series 
in 0" in which only the constant and hnear terms may depend upon {x,y), 
i.e., In0 = ?/,t)(7 + J2i=2 f^ii^)'^^ ~ ^^Z{x,t)- In such a case the term 
containing O in ( l66l) vanishes since by construction (O) = (o"0) = 0, confer 
(19), and one gets that 

^-^ = ^jdajdxdy'-m'>0. (67) 

It is for this reason that the only mixing measure considered in the body of 
this paper is the 0-th order So{t). Non-extensive quantities like Sr(cr, t) with 
r 7^ have been considered as possible generalizations of the Boltzmann 
entropy, see e.g. [28] and in the context of 2D fluid motion[2j , confer also 
reference [1] in which some useful comments are given. 



B 

In order to determine the family of probability densities (f)LB{cr,x,y,t) that 
reach the lower bound in i^S), i.e., (|Vln0LB|^) = iVc^l^ , we write 

ln0 =: ^bnix,y,t)a"', 

n=0 

and notice that, due to the normalization condition (11 II) . one has that 
exp {-bo{x, y,t)) = Jda exp f ^ 6„(a;, y, t)a'' 



It follows then that, 

V6o(x, y,t) = -J2 (^") V&n(a;, y, t), 

n=l 

and consequently that 

Vln0 = ^[a"-(O]V6„(x,i/,t), 

n=l 

and Vuj = J2 ((^ - ^) K - (^")]) ^bn{x, y, t). 



n=l 



It is now easy to see that if for n > 2 one has V6„ = then <^|Vln0| ) = 
((cr — cj)^) |Vfei|^ and also Vcj = ((a — cj)^) V6i(x, y, t), i.e., iVujf = |V6i|^ 
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Therefore, in this special case with Vbn = for n > 2 the equahty |Vci;| = 
K2 (|Vln0|^) holds. We have shown then that (f)LB{<^,x,y,t) can be identi- 
fied with (f)*{a,x,y,t) as given in Equation ( l3T|) with 

X{x,y,t) := bi{x,y,t), 
and := ^yU„(t)cr'", V/i„ = 0, 

n=2 

and that it attains the lower bound in ( 14^ . i.e., 



(|Vln0iB|'> 



|Vcu| 



C First and second variation of Jdxdy ( | V In 



In this Appendix we investigate the extrema of jdxdy (|Vln(/)|^) when 
satisfies ffTTl). ffT2l) as well as 



ci^(t){(J^ X, y, t) = m2{x, y, t). 



The quantity we have to vary, let us call it T^cj)), is then 



1 ^ 
T(0) :=-|V0|V^A„(x,2/,t)a^ 



n=0 



where Xn{x,y,t) are the Langrange multipliers corresponding to the local- 
moments constraints mQ{x,y,t) = 1, mi{x,y,t) = uj{x,y,t) and m2{x,y,t). 
Up to second-order terms in 6(f) we have 

r(0 + 50) - r(0) = 2V ■ (50V In 0) + 61T + 62T, 



with (5iT := 



and 62T :- 



2 

J_|V0|2-2^ + X^A„(x,i/K 



n=0 



V0 



The first term is a total divergence that, upon integration over {x,y), leads 
to a vanishing boundary term. The extrema are determined by 5iT = and 
since doT > it follows that all extrema are minima of Jdxdy (|Vln0|^) 
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under the given constraints. In order to solve 6iT = it is convenient to 
introduce e{a,x,y,t) := 0^^^ and noticing then that 

1 /A0_1^^ 

one sees that 5iT = is equivalent to 

Ae 



\ (^A„(a;,?/)(T" I e. 

\n=0 / 



4 

Let us denote the minimizer of T(0) by (t>ext{o', x, y, t), i.e., 6iT{(f)ext) = 0. As 
one can check, this minimizer can be written as 

(pexticr, X, y, t) = exp t) + bo{x, y, t) + y, t)a] , 

where exp [—bo{x, y, t)] := jda exp [/i(cr, t) + y, t)a] , 

which can thus be identified with (f)*{a,x,y,t) as given in Equation (IHTl) . In 
this extremal case, the Lagrange multipliers are, 

Ao = ^'|Vx|'-2a;Ax-2Vx- Vo; 

= [fsTa]"^ [{uJ^ - '2K2) I Vcj|^ - 27^2^ (Auj - Vcu ■ V In K2)] , 
Ai = 4Ax-2a;|Vx|' 

= [Ka]"^ [2K2 (Aiu - Vcj ■ VlnfsTa) - 2uj |Vcu|^] , 
A2 = |Vx|' = [i^2]"' |Vcu|% 
A„(x, y,t) = 0, for n = 3,4, .... 

D Cases associated with finite-order diff^erental 
equations 

The differential equations (1301) and fl3^ are of infinite order, however, in 
some cases they reduce to finite-order differential equations in a non-trivial 
way. For example, when is such that 

:r = — r(T, 68 
da QhAqa) 

where Pn{x) and Qm{x) are polynomials in x of order iV and M, respec- 
tively, lim^^^o = lim^^^o QmI^;) and g is a parameter such that qa is 
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dimensionless. Then equation (1481) reads 



and using the recursion operator L defined in ( l33l) . it follows that 

2q^Pj,{qC)n{x) = QM{qC)x. 

The remarkable disappearance of the a— average can be traced back to the 
special x~dependence of the extremal distributions, confer equation (|33l) . 
Making now use of equation (15T1) . the linear differential equation (|52l) that 
Z{x) must satisfy when (J68ll holds follows, namely 

2q'P^{q^)^ = QM{q^) {xZ) • (69) 
dX dx dx 

In Subsection 186.31 with N = 2 and M = was considered. An 

example with Pn{x) = 1 is provided by a /i(cr) of the following form. 



In [1 - d^gV^] for < [rfg]- 
— oo for cr^ > [dq]^'^ , 



where d is a pure number. For finite d the vorticity distribution has a finite 
support, (T^ < while in the limit d ^ it approaches a Gaussian. 

Since 

d^d 2q'^cr 



2' 



da 1 — (Pq'^a 

the differential equation satisfied by the corresponding partition function 
is 



dx 

(PZa 2 ,^ ,_o, dZd Zd 



1 j2 2 " 



ie — 't^ I (11 d^"^) = 



This is a modified Bessel equation of fractional order. The corresponding 
boundary conditions are Z^iS^) = 1 and Z'^{0) = 0. The solution is: 

\xj qd 
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where Ya{x/<id) is a modified Bessel function of fractional order and a : = 
{(f + 2)/2d?. Indeed, since in the present case exp //^(cr) = [1 — (fq^cr'^]^^'^ 
for (T^ < 1/ {q^dP) and otherwise, we have that 

/l/^f^ 2 
■l/gd 

corresponding to the integral representation of the modified Bessel functions. 
As one can check lim^_+oo ^(x) = [q'c?]~^and lim^^o '^^/'^X = l/g(3(i^+2).The 
graph in Fig. 2 shows the corresponding f2(x) for the cases d = 1 and o? = 3. 
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